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Abstract 



The theory of orbital magnetism in disordered metals is reviewed, and 
extended to include a broad range of temperatures and fields. Sample-to- 
sample fluctuations in the orbital magnetic susceptibility are studied. In 
a given sample these fluctuations manifest themselves in aperiodic sample- 
specific oscillations of susceptibility and magnetization, when the strength of 
the magnetic field is changed. 
1. Introduction 

A degenerate electron gas in a macroscopic metal exhibits weak orbital 
magnetism (the Landau diamagnetism) [1]. There are two factors that can 
significantly enhance the orbital magnetic response: finite system size and 
electron-electron interactions. This paper is devoted primarily to reviewing 
and extending recent work on the finite size (i.e. mesoscopic) effects on 
orbital magnetism of disordered metals (some of the results obtained prior 
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to 1993 have been reviewed in [2]). The effect of interactions will be briefly 
addressed at the end of the paper. 

In a weakly disordered metal the elastic mean free path I is much larger 
than the inverse Fermi number hp 1 but much smaller than the sample size 
L. Under these conditions electrons propagate by a diffusion process, with 
a diffusion coefficient D = VF^/d, where vf is the Fermi velocity and d is 
the sample dimensionality. In addition to the three length scales (hp 1 , £, L) 
which characterize the sample, there are two more lengths, L T = (hD/T)% 
and L B = (hc/eB)^, which specify the temperature T and the external 
magnetic field B. We assume that the inelastic scattering (or the phase 
breaking) length is larger than Lp and, thus, irrelevant for our purpose. 

Whether a sample of a given size L can be considered as macroscopic, 
depends on the temperature. If T is smaller than the Thouless energy E c ~ 
hD/L 2 (i.e. Lp > L), the sample forms a single coherent unit and can exhibit 
large mesoscopic effects. The orbital magnetic susceptibility x fluctuates 
from sample to sample and the typical fluctuation, (A^ 2 )^, can exceed the 
average value (x) by a large factor (angular brackets denote averaging over 
an ensemble of macroscopically identical samples). Mesoscopic effects in 
orbital magnetism are not restricted to T < E c but persist to much higher 
temperatures, although their magnitude gradually decreases [3-6]. A detailed 
study of (A% 2 )5, in the temperature regime ^CTc Hvf/£, was made in 
[5,6]. Some of these results will be re-derived below. The present treatment, 
however, allows us to consider a much broader range of temperatures and 
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fields. In particular, temperatures larger than Hvf/£ (i.e. Lt < £) and fields 
stronger than hc£ 2 /e (i.e. Lb < (!) will be discussed. 

In Sec. 2 we explain the phenomenon of mesoscopic orbital magnetism 
using an example of a clean chaotic cavity. A clear understanding of this 
simple system is helpful for the qualitative considerations of Sec. 3. A quan- 
titative theory for (Ax 2 )^ of a disordered metal is set up in Sec. 4, and its 
consequences are examined in Sec. 5, for various temperatures and fields. 
Sec. 6 is devoted to a brief discussion of the interaction effects. Conclusions 
are summarized in Sec. 7. 



2. Chaotic cavities 

Before turning to disordered metals, let us briefly discuss the simpler case 
of a clean chaotic cavity. We follow below the presentation in [2]. (Similar 
considerations for persistent currents appear in [7,8].) The two relevant ener- 
gies are the level spacing A and the "inverse time of flight" across the cavity, 
T ~ Hvf/L. The mesoscopic regime corresponds to the temperature range 
A <C T<X. In this regime the thermodynamics of the system is determined 
by the density of states pr(e), smoothed over an interval of order T near 
energy e. The main point is that pr contains a small oscillatory component, 
5pr, which is due to the shortest periodic orbit [9] 

1 (2ixBL 2 \ 
5p r (e,B) ~ -cos(A;L + 7) cos ( — J . (1) 

where k = h~ 1 \/2me, O = 2nhc/e, and the magnetic field B is assumed to 
be weak, so that the cyclotron radius R c = VFfnc/eB is larger than the size 



3 



L of the cavity. The oscillatory term in the density of states contributes a 
small oscillatory correction 5Q to the thermodynamic potential: 



which, after integrating by parts twice and keeping the leading term in the 
large parameter kpL, gives 



It is this correction that, due to its sensitivity to the magnetic field, dominates 
the differential magnetic susceptibility 



where xo — e 2 kfr 2 /mc 2 is the absolute value of the Landau susceptibility, in 
d dimensions. 

Thus, the orbital susceptibility of a two-dimensional electron gas confined 
to a chaotic cavity (a quantum dot) is enhanced by a factor fcpL, as compared 
to Xo- Detailed semiclassical computations of the orbital magnetic suscepti- 
bility for a generic chaotic system have been performed by Agam [10], and for 
a "stadium billiard" by Raveh [6]. Similar semiclassical considerations apply 
also to integrable systems [11-13] where, due to existence of families of peri- 
odic orbits, the enhancement is even larger. In all clean systems mesoscopic 
effects rapidly disappear when temperature is raised above T, i.e. for T ^> T 
the system becomes macroscopic and its orbital susceptibility assumes the 
Landau value. This is in contrast to the situation in disordered metals, where 
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the effect decays with temperature rather gradually. 

3. Disordered Metals: Qualitative Considerations 

The following qualitative discussion is confined to two-dimensional sam- 
ples. In a disordered metal electrons propagate diffusively, so that the en- 
ergy T of the previous section should be replaced by the Thouless energy 
E c = hD/L 2 . For T<E C , a change 5B ~ O /L 2 in a magnetic field produces 
a change 5Q ~ E c in the thermodynamic potential. The typical magnetic 
susceptibility Xt yP of a given sample can have either sign, and its magnitude 
is of order 5Vt/L 2 5B 2 ~ Xo^f^- 

The large sample-to-sample fluctuations cancel out, to a large extent, if 
one averages over many samples, with different impurity arrangements. The 
resulting average value (%) is much smaller than x typ and depends on whether 
the samples are in contact with particle reservoir (grand canonical system) or 
are isolated (a canonical one). In the former case (x) ~ — Xo = — e 2 /127rmc 2 , 
whereas in the latter there is an additional paramagnetic contribution which, 
for T <C E c , is the dominant one [4,14]. We will concentrate on the typical 
susceptibility, or, more precisely, on its fluctuating part (Ax 2 )^ = $X- This 
quantity is insensitive to the type of thermodynamic ensemble. 

Thus, for T^E C , 5x — Xo&f^- Let us see how 5x changes when T is raised 
above E c . For E c < T < (hv F /£) = h/r (i.e. i < L T < L), one can view 
the sample as made up of (L/L T ) 2 boxes, of size L T each. The above estimate 
is applicable for each box separately, i.e. Sxbox — Xo^f^ an d the fluctuating 
part of the magnetic moment of a box is SM Box ~ x<^f^BL\. Summing 
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over all boxes, with random signs, one obtains the typical fluctuation of 
the magnetic moment of the entire sample, SM ~ xokp£BL T (L/ L T ), i.e. 
5x/xo — k F £(L T /L) ~ T~5. This ratio is smaller than k F £ but can be 
larger than unity. Even if this relative fluctuation is small, it can still be of 
interest because of its non-trivial temperature and field dependence (strong 
field effects will be discussed later). 

The picture changes, if the temperature is raised further and becomes 
larger than h/r (i.e. L T < £)■ The relevant trajectories, which dominate 
the magnetic response, are now of ballistic (rather than diffusive) nature, i.e. 
shorter than t. Only short trajectories, of an area S^hvp/T) 2 = £ T , can 
efficiently respond to the magnetic flux. Let us denote by Wt the probability 
for having this kind of a trajectory, in a "box" of size It- The probability 
P(S)dS for a closed ballistic trajectory, of an area between S and S + dS, 
was estimated in [15] as P(S)dS ~ £- 3 VSdS, so that W T ~ {£ T /(-f- The 
fluctuating part of the susceptibility 5x of the sample is estimated in a way 
similar to that used above for the diffusive regime: The sample is partitioned 
into (L/£ T ) 2 boxes; a box contains, with probability Wt, a ballistic trajectory 
of an area of order £ T which contributes a magnetic moment |5M Bo x| — 
Xo^f^t B£ T (the factor kp£T is the enhancement factor discussed in Sec. 2). 
Adding up the (random) contributions of all boxes, and remembering that 
only a fraction Wt of all boxes does carry magnetic moments, one obtains 
5M ~ \5M Box \(W T L 2 /£ T )^ i.e. 5 X ^ Xok F £(£/L)(£ T /£)h Thus, the T~\ 
dependence of the diffusion regime (E c < T < h/r) crosses over to a T~2 
dependence in the ballistic regime (T > h/r). 
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The arguments so far were restricted to weak magnetic fields, i.e. to the 
regime of linear response. The region of validity of linear response depends, of 
course, on temperature. For instance, in the diffusion regime linear response 
requires L B > L T , i.e. fiuj c < T/k F £, where uj c is the cyclotron frequency. For 
stronger fields various non-linear effects in susceptibility and magnetization 
show up. In particular, for Lb < Lt, one may expect aperiodic, sample- 
specific fluctuations in magnetization [6]: every time when the magnetic 
field acquires an increment AB ~ 0o/^r, an effectively "new sample" is cre- 
ated and, thus, the differential susceptibility changes by an amount of order 
^Xo^f^{Lt/L). These sample-specific fluctuations resemble the universal 
conductance fluctuations (see [16] for an extensive discussion), although the 
relative fluctuation in susceptibility is much larger than in the conductance. 

4. Disordered Metals: Theory 

The grand potential, for a specific sample, is 



where V = L d is the sample volume, z/(e, B) is the density of states at field 
B, per unit volume, and 



The sample magnetization is M = —V^^Vl/dB)^^ and its differential 
susceptibility is x — dM/dB. The variance, (Ax 2 ), in the ensemble of 
random samples is therefore 




(5) 
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de^F^ti - n)F T (e 2 - jj)K{e l , e 2 ; B u B 2 ), (7) 
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where K(ei, e 2 ; B±, B 2 ) = (Ai/(ei, Bi)Ai/(e 2 , B 2 )) is the correlation function 
for the density of states, at two different energies and for two close values of 
the magnetic field. 

The density of states correlation function for diffusive metals was first con- 
sidered by Altshuler and Shklovskii [17]. In the context of orbital magnetism 
it has been studied by a number of authors [4-6,14,18]. The computations 
are based on writing the density of states in terms of retarded and advanced 
Green's functions: 

Z/(C ' B) = ~ S 2^V Tr [Gr ^ B) ~ Ga ^ B)] ' (8) 

where the factor s = 2 accounts for spin degeneracy. In what follows it is 
useful to use the approach of Altland and Gefen [18] and to write 



G R , A (e,B) = ^-\n(e±t V -H -V) = 



d_ 
d~e 



00 1 

1 / ~(0) 



He±ir}-H ) + J2-(G^ A V) 



n=i n 



(9) 



where H = (l/2m)(p — f^4) 2 is the unperturbed Hamiltonian, G$ A are the 
corresponding Green's functions and V is the random potential. Eq.(9) is 
written in an operator form, rather than in the coordinate representation. 
Averaging of a product Gr(6i, Bi)Ga(^2, B 2 ) results in the usual "diffuson" 
and "cooperon" series. The contribution of the n-th term of the cooperon 
(diffuson) series to the correlation function K(e±, e 2 ; Bi, B 2 ) is: 



K (n) = S 2 d 2 

± " 2tt 2 V 2 de x de 2 



-ReTrPl{e x -e 2 ,B ± ) 



(10) 
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where sign +(— ) corresponds to cooperon (diffuson), and B± — B\ ± B 2 . 
The operator P + , in coordinate representation, is written as 

P + (f, e l - 62, B + ) = JL-(G R (?, 7; €l , Bi))(Ga{t, 7; c 2 , B 2 )) , (11) 

where u is the free electron density of states (per spin), for an infinite system 
at B = 0. The expression for P-(r, f'; e\ — e 2 , B ) is obtained from Eq. (11) 
by interchanging r and f", in the advanced Green's function. For classically 
weak magnetic fields, i.e. when uj c t <C 1, 

(G RA (r, 6, B)} = (G RA (r, ?• e,B = 0)) exp ^ jf" A ■ dlj , (12) 

where the integration is along a straight line connecting r to f". 

To obtain the full function K(ei,e 2 ; B i} B 2 ) one has to sum K±^ over 
n. Since, eventually, derivatives with respect to Bi,B 2 are to be taken 
[Eq. (7)], one is interested only in the independent part of K. The (sin- 
gle) term with n — 1 and the two terms with n = 2 (one for diffuson, one 
for cooperon) do not depend on B. This happens because the independent 
phase factor drops out in a diagonal term (Gii(r,r)), as well as in the prod- 
uct (G R (f,f'))(G R (f',r)) (and similarly for the advanced Green's functions). 
Thus, the final expression for the -B-dependent part, AK, of K is: 

2 Q2 » 1 

AK (£1 , f2l Bl , B2) = ^^-Re g B + H 

+A;(e,-e 2 ,B_)) , (13) 

where A a (ei — e 2 , -B±) is the a-th eigenvalue of the operators P±. 

Eq. (13), supplemented by the definition (11) of the operators P±, forms 
the basis for a quantitative treatment of orbital magnetism in disordered 
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metals, in a broad range of temperatures and fields. The treatment is not 
restricted [18] to the standard diffusion approximation for cooperon and dif- 
fuson. The situation is the same as in the treatment of the weak localization 
correction for the conductivity. This correction is expressed in terms of the 
eigenvalues of the operator P + and, in order to obtain reasonably accurate 
quantitative results in a broad range of magnetic fields, one has to evaluate 
these eigenvalues beyond the diffusion approximation [15, 19-21]. 

Physical quantities of interest can be expressed in terms of AK, given in 
Eq. (13). Along with the variance 

(Ax 2 ) = \™ BdB T dB2 J f de 1 de 2 F T (e 1 - t i)F T {e 2 - f i)AK(e 1 ,e 2 ;B 1 ,B 2 ) ,(14) 

B 2 ^B 1 2 

one can consider various correlation functions, such as (Ax(B)Ax(B + AB)) 
or (Ax{T)Ax{T + AT)). The first function is defined by Eq. (14) with 
B 1 — > B, B 2 — > B + AB. To obtain the second function one has to replace 
one of the F T -factors in Eq. (14) by F T+AT . Other correlation functions, for 
susceptibility or magnetization, can be defined in a similar way. 

Integration in Eq. (14) can be carried out. After integration by parts, the 
derivatives contained in AK act on the F T -factors, giving 

^Me-n). (15) 

The subsequent integration reduces to the Matsubara sums over the poles 
of fx(e — //). This introduces a factor (2nT) 2 and changes the argument 
ei — e 2 of the functions under the integral to e p + e q where e p = mT{2p + 1), 
e q = mT(2q + 1) and p,q = 0,1, The final result is (the double sum 



d 
de 



F T {e-n) 



exp( 



T 



) + l 



10 



over fermionic frequencies can be replaced by a single sum over a bosonic 
frequency) : 

+\Z(e p + e q ,B_)\ (16) 

In the next Section we study (Ax 2 ), for various temperatures and fields, and 
briefly discuss some correlation functions. 

5. Results and Discussion 

One should distinguish between weak magnetic fields, L B > L T (linear 
response) and stronger fields, when the dependence of magnetization on B 
becomes non-linear. Note that Lb > L? is a sufficient condition for linear re- 
sponse, i.e. there is no requirement that Lb should be larger than the sample 
size L (recall that we are interested in the "high temperature" case, L T < L). 
In this respect the situation is like in the thermodynamic perturbation the- 
ory [1], where finite temperature stabilizes the perturbative treatment and 
enlarges its range of validity. Thus, while computing the eigenvalues of the 
operators P±, we can assume L B < L. Moreover, it is sufficient to find only 
the eigenvalues of P + . This is because the function P differs from P + only 
by its argument, containing (P>i — B 2 ) instead of (Bi + B 2 ). In particular, 
since in the definition of (Ax 2 ) the difference (Bi — B 2 ) is infinitesimal, the 
contribution of P_ to (Ax 2 ) is the same as of P + in the limit of small B. 

Eigenvalues of the operator P + are well known, in connection with the 
weak localization problem [20] (a closely related operator appears also in the 
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theory of superconductivity [22,23]). For the two-dimensional case, which is 
under consideration here, the N-th eigenvalue of P+(e\ — €2, B + ) is 



(17) 



IIjv(ei - e 2 , 5+) = — dxL N (x 2 ) exp 

€ JO 

and each eigenvalue is degenerate (L 2 /nL 2 + ) times. L N is the Laguerre poly- 
nomial (N — 0, 1, . . .) and = 2hc/eB + . Thus, the contribution of P + to 
the variance of susceptibility, Eq. (16), is: 



(AX 2 )+ = —T7- Jim 



L 4 Bi - bQ2 Bi Q2 B2 

B 2 -> B 



T 2 00 00 00 1 

7rlj + p,g=0 n=3 7V=0 Ai 



(18) 



where 



nAf(e p + e^, B + ) = -j- dxL N (x 2 ) exp I — — - 



-x 



, 27rTr , 1 \ 



To obtain the contribution (Ax 2 )~, corresponding to P_, one has to replace 
£>+ by P_, in the argument of IIjv and in the definition of L + . Thus, for any 
value of the actual magnetic field B (satisfying, of course, u c t <C 1), (A% 2 )- 
is given by the weak field limit of (A% 2 )+- 

For Lb-iLt ^> t (the diffusion regime), the integral in (18) is dominated 
by small x, so that L N (x 2 ) ~ 1 — Nx 2 , and 

2£ 2 ,„ 1, 2vrTr 

which holds for N <C (L + /£) 2 and p,q <C U/Tt. Under these conditions 
the sum over n in Eq. (18) can be approximated by — ln(l — H N ) and the 
expression for the typical value (A% 2 )^ = 5x can be reduced to the one given 
in [5]. In the limit of weak field 

5 X =l-Q7Xo^fk F e , (21) 
12 



nSf) ^l--—(N+-)-^j-{ P + q + l) , (20) 



where xo = e 2 /127rmc 2 . 

It is known from the previous studies [15,18-21] that calculations based on 
the diffusion approximation for IT7V are quantitatively accurate only when the 
ratio £/L B is extremely small. In order to study larger fields or higher tem- 
peratures, in particular when Lb or become smaller than £ (the ballistic 
regime), one must treat eigenvalues 11^ more accurately. Without resorting 
to numerics, one can use the approximate expression [20] 



which interpolates between the two regimes. In the ballistic regime, as op- 
posed to the case of diffusion, II jy is small so that it is sufficient to keep in 
Eq. (18) only the term n — 3. Note that the present approach, while giving 
the correct asymptotic dependence of 5\ on B and T, is not capable of pro- 
ducing the correct numerical coefficient in the ballistic regime. The reason is 
the that the ballistic regime is dominated by three-impurity diagrams, and 
there are several such diagrams besides the n = 3 diagram of the diffuson or 
cooperon series. 

Eq. (18), together with the approximate expression (22) for ITtv, enables 
one to consider a broad range of temperatures and fields. Let us first discuss 
the case of linear response, i.e. Lb > L?- In this explained above, 

(A% 2 )+ = (A% 2 )_ and does not depend on B. At low temperatures, L T 3> I, 
one is back to the diffusion regime, Eq. (21). At high temperatures (Lt < £, 
i.e. T > h/r) one enters the ballistic regime, where II jv is small. Keeping 
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(22) 
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only the term n = 3 in Eq. (18), one has to compute the sum 



E 



Ab{N + \) + t 2 



(23) 



N=0 

where b = (£/L + ) 2 and t = (2-kTt /%){p + q + 1). Although the sum can 
be computed by the Euler-McLaurin formula, it is instructive to extract 
the T-dependence of (Ax 2 ) by simple power counting. Since there are four 
derivatives with respect to B 1: B 2 in Eq. (18), one has to look for terms 
proportional to B+. One power of B + comes from the term L+ 2 , and the 
other three from the function defined by the sum in Eq. (23). Since this 
function is of the form t~ 3 f(b/t 2 ), it is clear that the term proportional to b 3 
will be multiplied by t~ 9 . The sum over p and q rapidly converges, so that 



(A X 2 ) * 



t 2 d A 

L 2 dB 4 



(24) 



i.e. 5x — X^k- F i{JL / 'L)(£ T / 'tj^ , in agreement with the arguments of Sec. 3. 

Consider now the opposite case, L B < L T) when (A% 2 )+ (but not (Ax 2 )-) 
begins to depend on B. Since (Ax 2 )+ decreases with B, it follows that for 
sufficiently strong fields (Ax 2 ) will reach half of its low-field value. In the 
diffusion regime, £ « L B « L T , (Ax 2 )+ decays rather slowly, on a typical 
scale of AB ~ ^jL\. For stronger fields, L B « £ « L T) the ballistic 
regime is reached. In this regime (Ax 2 )+ becomes negligible, compared to 
(Ax 2 )-, and keeps decreasing as L 7 B . Indeed, in this case the sum (23) is 
of order 6~ 3//2 , as long as p, q < Vb(h/27rTr), and the expression in the 
square brackets in Eq. (18) is estimated as [L 2 / ' L 2 + )b- z/2 b{h/2-KTr) 2 ~ y/B^. 
After differentiating four times with respect to B + , one obtains (Ax 2 )+ — 
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Xo{kF^) 2 {L 7 B /L 2 £ 5 ). Thus, just on the border of the ballistic regime (Lb — t), 
(Ax 2 }+ — Xo(^f^) 2 (^/L) 2 , i.e. is reduced by a factor (i/L T ) 2 as compared 
to its low-field value. 

Various correlation functions can be studied in a similar way. Consider 
the functions ((Ax(B)Ax(B + AB)) = C(B,AB). This function, as was 
mentioned in the previous Section, differs from (Ax 2 ) only by the limiting 
value of B 2 , which is now equal to B + AB. Therefore, the "cooperon part", 
C + , is a function of 2B + AB, whereas the "diffuson part", C_, is a func- 
tion of AB. Both C + and C_ decay with AB. In addition, C + decays as 
a function of B and becomes negligible for Lb <C Lt- Thus, for a fixed 
temperature (Lt £) there should exist aperiodic oscillations of magnetic 
susceptibility, in an individual (typical) sample, with a characteristic period 
AB ~ 4>q/L t and amplitude 5x — Xq^f^(L t /L). Since C_(AB) does not 
depend on B, these oscillations will persist into the ballistic regime, Lb <C £, 
and are limited only by the condition oj c t <^ 1. The possible existence of such 
oscillations, in the diffusion regime, was already mentioned in [6], although 
it was not realized that these oscillations do not decay even in the ballistic 
regime. Similar oscillations should exist, when the temperature or the chem- 
ical potential are changed by an amount of order T (or E c , if T < E c ). Such 
a change produces a "new sample" , and, thus, a corresponding change in the 
susceptibility. 

6. Interactions 
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Interactions can have a profound effect on orbital magnetism. Long ago 
Aslamazov and Larkin [24] pointed out the existence of an interaction induced 
paramagnetism. This phenomenon is similar to the interaction induced dia- 
magnetism in superconductors above the critical temperature. Since, how- 
ever, in a normal metal the interaction is repulsive, one ends up with a 
paramagnetic term. In two dimensions this term supercedes Landau dia- 
magnetism of non-interacting electrons. 

Static disorder modifies but does not destroy the effect. The combined 
effect of interactions and disorder was studied in detail in [25,26], where the 
average susceptibility for an ensemble of disordered samples was calculated 
(for a semiclassical derivation see [27]). The effect is due to the interaction (in 
the presence of disorder) between a pair of electrons in the Cooper channel. 
This interaction, at zero magnetic field, produces a small correction to the 
single particle density of states. In two dimensions the correction, at the 
Fermi energy /i, is: 

1 , ln(T /T) 

where T = /zexp(l/A ) and A is the effective electron-electron interac- 
tion constant. If A were much smaller than l/ln(/x/T), one could expand 
the ratio of the logarithms in Ao, to obtain the first order result bv c ~ 
(Xo/hD) XniTr/h). Since, however, the aforementioned condition is violated, 
one should use the renormalized interaction in the Cooper channel which 
results in the double logarithm in Eq. (25), with T roughly equal to //. 

The magnitude of the interaction induced paramagnetic susceptibility (in 
the linear response regime) can be estimated by using the general relation 
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between susceptibility \ and the value of the effective Bohr magneton (3 e & for 
particles in question, namely: x — Pls v i where v is the corresponding density 
of states. The "particles" (cooperons) have charge 2e and "mass" h/2D, so 
that P eS ~ eD/c (larger by a factor k F £ than the electron Bohr magneton 
eh/2mc). The corresponding density of states is the correction 5v c , Eq. (25). 
Thus, the average paramagnetic susceptibility is 

W^^^Si- (26) 

Note that a similar estimate could have been made for the typical sus- 
ceptibility of a mesoscopic sample, in the absence of interaction (Sec. 3). In- 
deed, multiplying P^ s by the mesoscopic correction to the density of states, 
\8v\ ~ 1/hD, gives x — kF^Xo- Let us emphasize that there is a big difference 
between the mesoscopic correction 5u and the interaction induced correction 
5v c . The former is a finite size effect, sensitive to temperature and rapidly 
oscillating with /i. The latter is a robust, intrinsic property of the interact- 
ing system, so that the interaction induced paramagnetism is not at all a 
mesoscopic effect. 

The large paramagnetic susceptibility remains, of course, present also in 
the mesoscopic samples considered in Sec. 3-5. For T > E c the average sus- 
ceptibility (x) is already close to its macroscopic limit, Eq. (26). (For T <^ E c 
some finite size effects in (x) show up [4].) Thus, the mesoscopic fluctuations 
discussed in the present paper will occur on the background of (x P ). These 
fluctuations do not seem to be affected by interactions in a significant way, 
similarly to the case of the universal conductance fluctuations, as discussed 
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in [16]. However, the full picture of the mesoscopic magnetism in disordered 
metals, in the presence of interactions and for arbitrary temperatures and 
fields, has not yet emerged. 

7. Conclusions 

The theory of mesoscopic effects in orbital magnetism of disordered met- 
als was reviewed and extended to higher temperatures and fields. The effects 
decay with temperature rather gradually and persist to temperatures much 
higher than the Thouless energy E c . One manifestation of mesoscopic orbital 
magnetism is aperiodic, sample-specific oscillations in magnetization and sus- 
ceptibility under changes of various external factors - temperature, chemical 
potential or external magnetic field. The oscillations persist to quite strong 
fields, limited only by the condition oo c t <C I. The amplitude of oscillations 
in the susceptibility is of order Xo^f^ at T ~ E c , and slowly decreases at 
higher temperatures. 

One should distinguish between the average susceptibility, measured on 
a large ensemble of macroscopically identical samples, and the typical sus- 
ceptibility measured on an individual sample. I am not aware of any mea- 
surements of orbital magnetism in disordered mesoscopic metals. The sim- 
ilar phenomenon of orbital magnetic response to an Aharonov-Bohm flux 
has been investigated, to some extent, experimentally [28-31]. In particular, 
measurements on single rings, disordered [29] and quasiballistic [30], were 
made. Since, at low temperatures and weak fields, the magnetic moment of 
a ring is about the same as of a singly-connected sample of a comparable size, 
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the mesoscopic fluctuations studied in the present paper should be amenable 
to observation. 

Some aspects of the phenomenon of mesoscopic orbital magnetism were 
not covered in this work. We have not discussed mesoscopic effects for fields 
uj c t > 1, i.e. in the regime of the de-Haas- van Alphen oscillations or in 
the quantum Hall regime. Another omission is the mesoscopic magnetism 
in three dimensions, where a qualitatively new effect shows up. The sus- 
ceptibility of an individual sample becomes a tensor Xa/3 (ct,/3 = x,y,z), 
with non-zero off diagonal components, as opposed to the ensemble averaged 
susceptibility. This happens because only after averaging over the impuri- 
ties the sample becomes spatially homogeneous. In an individual sample the 
direction of magnetization does not coincide with the direction of the mag- 
netic field. For instance, for a sample of a spherical shape, the ratio between 
(Ax^ 2 )a and (Axj, z )^, is l/VS [5,6]. It would be of interest to investigate 
different sample shapes and to extend the treatment beyond linear response. 
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